SINP/TNP/2008/16 



Bound States in Gapped Graphene with Impurities : Effective Low-Energy 
Description of Short-Range Interactions 

Kumar S. Gupta* 

"qq" ■ Theory Division, Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Calcutta - 700064, India 

| Siddhartha Sen t 

■ School of Mathematical Sciences, VCD, Belfield, Dublin 4, Ireland and 

. Department of Theoretical Physics, Indian Association for the Cultivation of Science, Calcutta - 700032, India 

. Q , We obtain a novel bound state spectrum of the low energy excitations near the Fermi points of 

\ • gapped graphene in the presence of a charge impurity. The effects of possible short range inter- 

' actions induced by the impurity are modelled by suitable boundary conditions. The spectrum in 

the subcritical region of the effective Coulomb coupling is labelled by a parameter which character- 
izes the boundary conditions and determines the inequivalent quantizations of the system. In the 
supercritical region we obtain a renormalization group flow for the effective Coulomb coupling. 
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I. INTRODUCTION 



CO 

!> ' Graphene monolayers consist of carbon atoms arranged in a honeycomb lattice. In the neighbourhood of Fermi 
points, the low energy excitations in graphene can be described by a two dimensional massless Dirac equation The 
effects of charge impurities in graphene have to be analyzed separately in two regions, subcritical and supercritical, 
depending on the strength of the Coulomb interaction. In the supercritical region, where the effective Coulomb 
• ■ strength exceeds a certain critical value, the massless Dirac equation admits bound states [2], y, |J, [5j . In the subcritical 
region this does not happen, which is a manifestation of the Klein paradox. 

If the exact honeycomb lattice symmetry in graphene is partially broken, possibly due to the presence of an impurity, 
a Dirac mass for the excitations can be generated. The effect of a charge impurity in gapped graphene with massive 
Dirac excitations has been analyzed in 0, @, @, Q, where it was assumed that the impurity provides an axially 
symmetric Coulomb interaction. It is expected that such an impurity may also induce other short-range or singular 
interactions, such as a delta function type potential. We have neither any detailed knowledge of such interaction 
terms, nor is it practical to include them in the Dirac Hamiltonian, which is valid only in the long wavelength limit. 
We can however still model the combined effect of these additional short range interactions on the long wavelength 
dynamics through the choice of suitable boundary conditions [Icj . 

In this paper we shall analyze the effects of these boundary conditions on the spectrum of the massive Dirac equation 
in the presence of a charge impurity. This approach, where self-adjointness is taken as the guiding principle for deter- 
mining the allowed boundary conditions, has been shown by Jackiw to yield a reliable description singular potential 
such as a delta function For instance, application of this approach to singular short range interactions has led to 
the possibility of novel bound states in fermionic [H, [TH, 0, [TH, HU and anyonic systems [l?], [H[ , molecular physics 
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1C | and integrable models [19|, |20| . This technique is particularly relevant for systems having scaling interactions 



22, 23, 24], a property present in the screening effect of the Coulomb potential in graphen e 125, 1261 I27I. l28l. |2 



Moreover, this approach has already been used to study certain topological defects in graphene [3fj, l3l|. It is thus of 
interest to explore the physical effects of generalized boundary conditions in gapped graphene with a Coulomb charge 
impurity. 

It should be noted that the analysis of the bound states in 0, @] assumes that the wavefunctions vanish at the 
location of the impurity, whereas we only require square-integrabilty for the wave functions. The boundary conditions 
that are consistent with this requirement introduce an extra parameter, which has physical implications. In particular, 
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the result for the spectrum obtained here is in general different from that obtained in [6|, |7[ ■ The parameter which 
appears in the boundary condition characterizes the inequivalent sectors of the quantum theory. The appropriate 
theoretical description for graphene in the presence of a charge impurity can be settled experimentally possibly 
through the STM measurements of the local density of states (LDOS). 

It is generally believed that the long wavelength Dirac description is not applicable to graphene in the supercritical 
region, where the Dirac vacuum is expected to break down [7J. However certain features of the system may still 
be captured by the continuum description. For instance, numerical as well as semi-classical analysis Q of the 
supercritical region in the massless case exhibits a large number of bound states in graphene. An analytical prediction 
of these bound states can be obtained from the continuum massless Dirac description Q . Here we shall apply the 
Dirac picture to the supercritical region with a cutoff comparable to the lattice spacing. We propose a renormalization 
group analysis by keeping the observables fixed as a function of the cutoff (ill . [32| , to evaluate the corresponding /?- 
function. An alternative regularization scheme for analyzing supercritical charge impurities in graphene with a mass 
gap has been discussed in [8| . It is not necessary to introduce such a cutoff in the subcritical region as the excitation 
energies are low compared to the inverse lattice scales. 

This paper is organized as follows. In Section 2, we set up the Dirac equation for the problem. In Section 3 we 
discuss the generalized boundary conditions that follow from self-adjointness. In Section 4 we find the spectrum 
with these boundary conditions. In Section 5 we discuss the spectrum in the supercritical region and the associated 
renormalization group flow. We conclude the paper in Section 6. 



II. DIRAC EQUATION 

We start by considering the massive Dirac equation in a gapped graphene monolayer in the presence of a charge 
impurity and use the same conventions as in Novikov Q. The Dirac operator can be written as 

H = -i{a 1 d x + a 2 d y ) + ma 3 + V(r), (1) 

where r is the radial coordinate on the two dimensional x — y plane, er,;, i — 1,2,3 are the Pauli matrices, and m 
denotes the Dirac mass. We have chosen units such that the Fermi velocity v = 1 and the Planck's constant U = 1. 
Using these conventions, the Coulomb potential V(r) is given by 

V(r) = --, (2) 
r 

where we choose the impurity strength a > 0, signifying an attractive potential (25j. In addition, we assume that 
the effect of the charge impurity is such that it induces short range and possibly singular potentials, such as a delta 
function, whose detailed nature is not relevant. In our approach, we assume that the combined effect of these short 
range and possibly singular potentials can be modelled by imposing suitable boundary conditions on the wave function. 
The Dirac operator (JTJ) satisfies the eigenvalue equation 

= E^> (3) 

where E is the eigenvalue and 

Here ipi(r) and ipi(r) denote the radial part of the wavefunction and <f> denotes the angle in the x — y plane. 

In this paper we focus on the bound states of the Dirac equation {SJ, which satisfy \E\ < m. Consider the ansatz 

tpl(p) = Vm + Ee-^p^-^xiip), (5) 
i/} 2 (p) = Vm - Ee-%p"-*X2(p), (6) 

where p = 2jr, 7 = \/m 2 — E 2 , v — y/ j 2 — a 2 and j = k+ h. In this Section we shall deal with the subcritical region 
of the Coulomb potential, which is given by a < j for any j. Since the lowest value of j — ^, in the subcritical region 
the effective Coulomb strength must satisfy a < i. Furthermore, in terms of the variables P, Q defined by 



xi=P + Q, X2 = P-Q 



(7) 
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we get the equations 

M5)-(^%%2 + " + *)(5)-* 

which Hp defined above denotes the radial Dirac operator. These set of equations can be decoupled to give 

d 2 P , ,dP ( aE\ 
P^TT + (l + 2v-p)-r--(v P = (9) 



dp 2 dp \ 7 

d 2 Q ,„ „ ,dQ / aE 
dp 2 dp \ 7 



Thus we see that the functions P and Q satisfy the confluent hypergeometric equation [33j]. This equation has two 
linearly independent solutions, one of which is regular at the origin (denoted by M) while the other is regular at 
infinity (denoted by U). In Q, the boundary conditions were so chosen that the solutions were regular at the origin, 
which led to the wavefunctions 

. . . , — p ,, i r . ra«, . aE . , aE . , aE .. . 

i>x{p) = Vm + Ee-2p»-2[(j + )M(v ,l + 2v,p) + {v )M(l + v ,l + 2i/,p) (11) 

^7 7 7 

, , x / ^ p ,, i ,, mat. , aE . , aE . , aE 

Mp) = V^Ee-ip»-i[(j + )M{y ,l + 2v,p)-(v )M{l + v ,l + 2i/,p). (12) 

v 7 7 7 

The corresponding bound state spectrum was obtained in @, 0] as 

E pd = , (13) 

1 j- a 

with p = 0, 1, 2, for j > and p — 1,2,3,..., for j < 0. In the next Section we shall see that more general boundary 
conditions are possible which are consistent with all the requirements of quantum mechanics, leading to a different 
spectrum for the same Dirac operator. 



III. GENERALIZED BOUNDARY CONDITIONS 



In the usual description of quantum mechanics, it is assumed that the Hamiltonian is self-adjoint [34j |. so that the 
time evolution is unitary and the probabilities are conserved. In addition, for the bound states, the solutions should 
be square-integrable. In our search for the generalized boundary conditions, we shall be guided by these principles as 
formulated by von Neumann [34[. 

The Dirac operator H in Q consists of a radial and an angular part. The domain Y(cf>) on which the angular 
part of H acts is spanned by the periodic functions <&&(</>) , k G Z in (4). In what follows, we shall leave the angular 
wavefunctions and the corresponding boundary conditions unchanged. 

The radial part of the Dirac operator H is given by H p in ijHJl. It is symmetric (or Hermitian) in the domain 
Dq(H p ) = C§°(R + ) consisting of infinitely differentiable functions of compact support in the half line R + . The 
corresponding adjoint operator is denoted by H^, which, as a differential operator, has the same expression as H p in 
©, although its domain could be different. 

The domain Dq(H) of the full Dirac operator H is therefore given by Do(H) = Cfi°(R + ) ® Y (</)). Its adjoint 
operator has the same differential expression as H, although its domain could be different as well. Following 
von Neumann's approach [341 ] . in order to determine whether the full Dirac operator H is self-adjoint in its domain 
Dq(H), we consider the equation 

fl' t *± = ±i*±. (14) 

Let n+(n_) be the total number of square-integrable, linearly independent solutions of (|14p with the upper (lower) 
sign in the right hand side. The quantities n± are called the deficiency indices of H. In order to determine n+(n_), 
we consider the radial equation (jSJ) with E replaced by i). This will give the deficiency indices n± for H p . In 
terms of n± , H p can be classified as follows [3J| : 

1) H p is (essentially) self-adjoint in D (H p ) iff (n + ,n_) = (0,0). 

2) Hp is not self-adjoint in Dq(H p ) but admits self-adjoint extensions iff n + = n_ = n(say) 0. 
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3) H p has no self-adjoint extensions if n + ^ ra_. 

In order to find the deficiency indices n± for H p , we need to solve for P± and Q± from the equations 

(15) 



<fP± , ,dP± 



dp 2 



dp 



oa . 
v }P± 

7±, 



d 2 Q d 
1 dp 2 



dp 



0. 



(16) 



which are obtained from (9) and (10) with E replaced everywhere with ±z and where j± = v M 2 + 1. The solutions 
we seek are such that when we reconstruct 



and subsequently obtain 



Xi± =P±+ Q±, X2± =P±-Q± 

ip 2 ± = VmT ie~% p"~^X2±(p), 

the functions %p\± and ip2± would be square integrable on R + with a measure pdp. 

We now proceed to find n + . In this case, a possible set of solutions of (15) and (16) are given by 

P+ = U\v- ™,l + 2v,p 

Q+ = Un. + v-™,l + 2v,p) . 

where U denotes a confluent hypergeometric function [33]. As p — > oo, 

-v+l~ 



P 



(17) 

(18) 
(19) 

(20) 
(21) 



(22) 
(23) 



Using (17), (18), (19), (22) and (23) we find that as p — > oo, V'l+ji/W — ► 0. Hence the functions ipi+,ip2+ are 
square integrable at infinity. 

Let us now consider the behaviour of these functions as p — > 0. For this we shall use the formula [33[ 



U(a,b,z) 



7T 



sin irb 



M(a,b,z) 



x _ b M{l + a-b,2-b,z) 



T(l + a-b)T{b) 



r(a)r(2-6) 



where as p — > 0, M(a, b, z) — ► 1. Using (20) and ([24]) we see that as p — ► 0, 



Q+ 



where a 



sin 7r(l+2i^) 



and 



A+ = 
C + = 



Y{-v-^-)T{l + 2v) 



a{A + - B + p~ 2v ), 
a(C+ - D + p~ 2v ), 



Y{v-^)Y{\-2v) 



1 



Y(l-u-^)T{l + 2u) 



D + = 



1 



T{l + v-^)T{l-2v) 



are constants depending on the system parameters. From the above relations we find that as p 

\il>i+\ 2 pdp — » j {c lP 2v +c 2 + c 3 p- 2 »)dp 
\*lj2+\ 2 pdp — > / (d lP 2v + d 2 + d 3 p- 2 »)dp, 



(24) 



(25) 
(26) 



(27) 
(28) 

(29) 
(30) 
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where Ci,d iy i = 1,2,3 are constants whose explicit forms are not relevant. Recall that v = yj j 2 — a 2 and that in 
the subcritical region, a < j. Hence v is a real positive quantity in the subcritical region. Then, from (29) and (30) 
we find that 0i+, 02+ are square integrable at the origin provided v < j. Thus we arrive at the conclusion that the 
functions 0i+, 02+ are square integrable everywhere provided < v < i. Alternately we can say that the deficiency 
index n+ = 1 when < v < |. A similar analysis shows that for this same range of u, n_ = 1 as well. We have thus 

shown that when < yj j 2 — a 2 < i , the massive Dirac operator for graphene in the subcritical region of the effective 
Coulomb coupling is not self-adjoint in Dq(H p ), but admits a one parameter family of self-adjoint extensions. 



IV. INEQUIVALENT SPECTRA 

We would now like to find the spectrum of the system in the range of j and the effective subcritical Coulomb 
strength a such that < v = yj 'j 2 — a 2 < i, where the Dirac operator admits a one-parameter family of self-adjoint 
extensions. The deficiency subspaces for the radial Dirac operator H p are spanned by the elements 



,±_ ifeJlVS?i^(^-Q±)J' 

The domain in which the Dirac operator is self-adjoint is then given by D z (H p ) = D (H p ) © {c(e^?? + + e _ ^^_)} 
where c is an arbitrary complex number and z S R mod 2n [HI]. Thus we have a one parameter family of self-adjoint 
extensions, labeled by a real parameter z. For each choice of the parameter z, we have a domain of self-adjointness of 
the radial Dirac operator defined by D z (H p ). An arbitrary element rj z € D z (H p ) can be written as 



Viz \ _ ( (e" "01+ + e " 0i_) 
) V (e^02+ + e"^0 2 -) 



We note that as p — > 0, 



r)u \ f v / m + ie 2 Tp"-2(P + + Q + ) + y/ m - ie~^ p v ~i (P_ + Q_) 

r\iz J C V v^T^Ie^ p v -^-{P+ - Q+) + Vm + ie^i p v ~h (p_ - Q_) 



(32) 



(33) 



where P_ and Q- denote the complex conjugates of P+ and Q+ in (25) and (26) respectively. 

We now proceed to find the spectrum of the system when the boundary conditions are governed by the domain 
D Z {H P ). A solution of the physical eigenvalue problem can be written as 

^=N( v ^ e :t p zi {p J ( il) (34) 



y/m-E e-%p y -2(P - Q) 

where the functions P and Q satisfy (9) and (10) respectively and N denotes the normalization. Solutions of (9) and 
(10) that are square integrable at infinity are given by 

P = u(i>-^ ) l + 2v,p\ (35) 
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= U[l + v-—,l + 2v,p\. (36) 



Using (24), (33) and (34) we see that in the limit as p — > 0, 



P — » a{A-Bp- 2v ) (37) 
Q — ► a(C-Dp- 2u ) (38) 



where a = -. — 71 . „ N and 

sm 7r{ l+2v) 



1 



A r(-i/-^)r(i + 2i/) B v{u-^)T{i-2u) (39) 

C = =r D = =r . (40) 

r(i - v - ^s)r(i + 2v) r(i + i/-^)r(i-2i/) v ; 



0.90 0.92 0.94 0.96 0.98 

E 



FIG. 1: (Color online) A typical plot of f(E) in (43) with m = 1, j = § and a = 1.46. The horizontal line corresponds to the 
right hand side of (43). It can be shifted up or down by changing the self-adjoint extension parameter z. 

Hence, as p — > 0, 



V-^ aN { ^~E [{A - OfT-h -(B- D)p—i] ) (41) 

The physical solution ip in (41) must belong to the domain of self-adjointness given by D z (H p ). In fact, behaviour 
of the elements of the domain D z (H p ) determine the boundary conditions for the system. If %ji e D z (H p ), then 
as p — ► 0, the coefficients of r u ~^ and r~~ u ~^ in (33) and (41) must match. Comparing such terms and defining 
yjm + i(A + + C+) = £ie idl and y/m + i{B+ + D+) = &e 102 , we obtain 

7 2 V A + C 6cos(^ + f) 



1 + APJ B + D £ 2 cos(6> 2 + §) ' 



(42) 



Using (39), (40) and (42) we get 
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= - r ( i-2,)r(i + ,-^)(i-,-^) ^ cos(fli + f) 

1 i • r(i + 2i/)r(i- t /-fs)(i + i/-^) &cos(^ + f)- 1 J 



Eqn. (43) determines the spectrum in terms of the system parameters and the self-adjoint extension parameter z. 
Each choice of z corresponds to a boundary condition described by the domain D z (H p ) and leads to an inequivalent 
quantum theory. It may be noted that the theory itself cannot predict which choice of the self-adjoint extension 
parameter will be realized in a given system and this parameter must be determined empirically. Equation (43) in 
general cannot be solved analytically. However, for the special choice of z = z\ such that Qi + 4f = f > we have 

Ea 

v = -n, n=l,2,3,.... (44) 

7 

This leads to the spectrum (fl3|) obtained by Novikov [1, 0| for < v < | . For another special choice of z = Z2 such 
that 6t + f = f , we get 

Ea 

-v = -n, n = l,2,3,.... (45) 

7 

For a general choice of z, the spectrum can be obtained numerically, and example of which is shown in Fig. 1. It 
may be noted that for a general choice of z, the spectrum we obtain from (43) is very different from that in (13), 
which was obtained previously [6j, |7|. The corresponding bound state wavefunctions (41) are square- i nteg rable, but 
not necessarily regular at the origin. This feature appears in graphene with topological defects as well [30l [3l|. 

V. SUPERCRITICAL REGION 



The supercritical region is defined by the effective Coulomb strength a 2 > j 2 for any j. This implies that in the 
supercritical region a > h and v = \J j 2 — a 2 = ±ip where p € R. We now proceed to investigate the supercritical 
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coupling region for the massive Dirac equation. A study of the massive Dirac equation with a regularized Coulomb 
potential has been discussed in Q. We shall focus on the excitations satisfying E 2 < m 2 and introduce a cutoff in 
the radial direction set by the lattice spacing in graphene. The cutoff restricts our analysis to the region where the 
Dirac equation holds. The corresponding eigenvalue problem is solved with a hard-core boundary condition given by 

ip(p = Pa) =0, po = 2r 7, (46) 

where ip is the two component wavefunction in the supercritical region and ro provides a cutoff in the radial direction. 
In this case, the upper component ipi(p) m (5) has two linearly independent solutions given by 



£(p) = Vm + Ee-'p^M (ip,- ^-,l + 2ip,,pj 
C(p) = Vm + Ee'zp^'^M y-ip - 1 - 2i[i, pj 



(47) 



(48) 



The general solution which satisfies the boundary condition (46) can be written as 

Mp) = E(P)C(A>) - ap)Z(Po)}- (49) 



As p — > oo, we get that 



ipi(p) — > vm + Ee 2 p 



r(i + gjM) r(i-2^) 
r(^-f) CM " ru-f)^ 



(50) 



In order for the wave function to be square integrable, the quantity in the parenthesis on the rhs of (50) must vanish. 
This gives the condition 

r(i + 2 ¥ )r(- ¥ -f) = ^) 

r(i-2i/i)r(*>-^) c(po)' 1 j 

Eqn. (51) follows as an exact consequence of our analysis. In order to gain some physical insight, we shall now use 
several approximations. The results derived below are therefore valid only in a qualitative fashion. First we assume 
that as the cutoff ro approaches the lattice spacing, the hypergeometric function M in (47) and (48) can be replaced 
approximately by 1. Strictly speaking this is ture when the cutoff tends to zero, but it is a reasonable approximation 
in the long wavelength limit. Second, we assume that E 2 < m 2 . In other words, we shall trust our results only for 
energy scales below the Dirac mass. Using these assumptions in (51), we get 

Y(-iu- —) 

where S is the argument of T(l — 2ip). In order to proceed, consider the energy scale such that ^ << 1. In this 
case, the l.h.s. of (52) is approximately independent of E and depends only on the system parameter p. Denoting 
the argument of T(—ip) by 9, we get 

1 0-5-2 P * 

7p = 7T- e " > 53 ) 
2r 



where p G Z and 7 P = ym 2 — E 2 . We can satisfy the requirement of ^ << 1 by restricting p suitably. We now 

assume that p, through its dependence on the effective Coulomb coupling a, is a function of the cutoff rg. We keep 
E p or equivalently 7 P invariant as the cutoff is varied, which gives the function as 

P{p) = -r ^ ~ -p 2 . (54) 
dr Q 

We see that the coupling p admits an ultraviolet stable fixed point at p = or equivalently at a = j for the angular 
momentum channel j, to which the system is expected to flow (Til l32l]. In particular, a tends to its critical value \ 
for the angular momentum channel i = \- 
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VI. CONCLUSION 

In this paper we have used the freedom to choose generalized boundary conditions to model the effects of short 
range interactions introduced by impurities in gapped graphene. We used this approach to investigate the super and 
subcritical regions for the effective Coulomb charge. For the subcritical region we found that the generalized boundary 
conditions introduce a self-adjoint extension parameter z which labels the different inequivalent quantizations for 
< \J j 2 — a 2 < \. For a specific choice of z, the result of [0, 0] can be recovered. In general the spectrum obtained is 
different. Thus an experimental approach for determining the appropriate choice of the boundary conditions labelled 
by z is in principle possible. 

For the supercritical region, the analysis suggests a renormalization group flow a — > j for the j angular momentum 
channel, where j is half integer. In particular, for j = \, the effective Coulomb coupling tends to its critical value 
a — i. This conclusion is valid in a very restricted region where — << 

In this paper we have considered only bound states. A similar analysis for the scattering sector would be relevant. 
In addition, the analysis of self-adjointness in bilayer graphene with impurities [35| would also be interesting. 
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